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ABSTRACT 

Between February 24-25, 2009, the EIS spectrometer onboard the Hinode spacecraft per- 
formed special "sit & stare" observations above the South polar coronal hole continuously over 
more than 22 hours. Spectra were acquired with the 1" slit placed off-limb covering altitudes up 
to 0.48 Rq (3.34 X 10^ Mm) above the Sun surface, in order to study with EIS the non-thermal 
spectral line broadenings. Spectral lines such as Fe xii /II 86.88, Fe xii /1193.51, Fe xu /1195.12 
and Fe xm /1202.04 are observed with good statistics up to high altitudes and they have been 
analyzed in this study. Results show that the FWHM of Fe xu /1 195. 12 line increases up to 
- 0.14 Rq, then decreases higher up. EIS stray light has been estimated and removed. Derived 
electron density and non-thermal velocity profiles have been used to estimate the total energy 
flux transported by Alfven waves off-hmb in polar coronal hole up to - 0.4 Rq. The computed 
Alfven wave energy flux density f„ progressively decays with altitude from /„■ - 1.2 • 10^ erg 
cm"^ s~' at 0.03 R© down to ~ 8.5 • 10^ erg cm~^ s~^ at 0.4 Rq, with an average energy 
decay rate l^f„ltJi - -4.5 • 10"^ erg cm"^ s~^ Hence, this result suggests energy deposition by 
Alfven waves in a polar coronal hole, thus providing a significant source for coronal heating. 

Subject headings: Sun: corona; Sun: oscillations; Sun: UV radiation; line: profiles; waves 



1. Introduction 



Many different physical processes have been proposed to explain the heating and acceleration of plasma 
emanating from polar coronal holes and resulting in the fast solar wind streams. Basically two competing 
theories have been proposed: the first one envisages ubiquitous magnetic reconnections occurring at the 
base of the corona (the so-called nano-flare heating, originally proposed by Parker.. 198 3. .1988) . leading 
to a dissipation of magnetic stresses produced by fieldlines footpoints motions at the coronal base (also 
referred to as direct current heating). The second theory invoques the propagation and dissipation along 
magnetic field lines of diff'erent types of MHD waves (e.g. magnetoacoustic, Alfven, ion-cyclotron waves). 
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also referred to as altern ating current heating (se e e.g. lKlimchukll2006l. for a review on tiie coronal heating 
problem). More recently. iDe Pontieu et all (1201 ih have suggested that also a ubiquitous coronal mass supply 
in which chromospheric plasma (i.e. spicules) is accelerated upward into the corona, and they provide 
constraints on the coronal heating mechanisms. For this reason, over the last decades many attempts have 
been performed to provide evidence for detection of waves in the corona and to quantify the energy provided 
by wave dissipation. A complete review of wave detections in the corona is far from the purposes of this 
work: here we briefly review main works published on the detection of waves as a broadening of spectral 
line profiles observed in corona above the limb. This eff'ect, which has been observed with spectroscopic, 
spectropolarimetric and interferometric measurements, is often interpreted as a signature of Alfven waves 
propagating undamped in a density stratified corona, because these waves are expected to be responsible 
for unresolved plasma motions occurring also with a component along the line of sight (both towards and 
away from the observer), hence leading to a spectral line profile broadening. If this interpretation is correct, 
a detection of a decrease, instead of an increase, of the line width above a given altitude, would imply a 
damping of Alfven waves, hence energy deposition in the corona. 

One of th e first evidences of increasin g line profile widths with altitude was obtained during a total 
solar eclipse by IChandrasekhar et all (Il99lb with a Fabry-Perot (FP) interferometer: they found a rise of the 
Fe XIV /15303 (A) coronal line-width from 0.07 Rq (measured from the solar limb) to 0.18-0.2 Rp, and then 



a tend ency to fall off with higher radial distance. More recently, this result was confirmed by iMierla et al. 



teOOSh . who studied the variation in the Fe xiv line width in the inner corona (above 0. 1 R©) wit h LASCQ/C 1 



(1996- 1998) data, acquired with a FP interferometer used as a narrow-passband, tunable filter. iMierla et al. 



(|2008h found that the Fe xiv Une width was roughly constant or increased with height up to around 0.3 Rq, 
and then it decreased. Even if the interpretation of these line broadenings as signatures of Alfven waves 
is not universally acc epted, recent results demonstrate that the lower corona is permeated by these waves. 
Tomczyk et al.l (|2007n identified low-frequency (<5mHz) propagating Alfvenic motions in the solar corona 
from linear polarization images of the Fe xm A1075 line. In the chromosphere Al fvenic motions have been 
observed by De Pontieu et al. ( 2007a ) with Hinode/SOT, and De Pontieu et al.l ( 2007bi) have studied their 
relationship with chromospheric spicules observed at solar limb. Moreover, it has been demonstrated by 
Hinode/XRT data that coronal X-ray polar jets occurs at a larger rate than previously reported (~ 10 events 
per hour) and that these jets are a ssociated with smal l scale mass ejections propagating at velocities close to 
what expected for Alfven waves dCirtain et al.ll2007l) . 



A very large number of published work were based on spectroscopic observations. First studies date 
back to the Skylab data and have shown an increase of line width as a fun ction of height above t he limb (lines 



emitt e d by chromospheric o r transition region io ns with T< 2.210^ K: iDoschek et alJI 19771 : iNicolas et al 
I1977I : iMariska et alTll978h . bassler et all (ll99oh found an increase in the line width up to 0.2 Rq above 



the limb by analyzing several line profiles from the off-limb quiet corona. Studies of the equatorial coronal 
regions from SOHO data (SUMER and CDS) have confirm e d the broadening of spectral line width as a 
function of h eight above the limb (e.g. 



Wilhehn et al 



polar regions, iBanerjee et al.l (|l998|) : Doyle et al.l (|l999|) . analyzing Si viii spectra by SUM ER, have found 



20041 : iHarrison et al.ll2002l : IPovle et al.lll998h . In 



an increase in the non-thermal velocity up to 0.25 R©. This result was confirmed also by IWilhelm et al. 
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(120041) through Mg x Une observations provided by SUMER in a polar coronal hole. lO'Shea et alj (l2003h 
observed, from spectra obtained by CDS, a decrease of the Mg x line width at 0.15 Rq, where the dominant 
excitation changes from coUisional to radiative. 

While there is a general consensus that EUV line profiles broden in polar coronal holes moving ofF-limb 
up to ~ 0.15 - 0.20 Rq, the line width behaviour above that altitude is not well determined . For instance, 
some authors found that th e line width profile shows a "platea u" above a given al t itude (e.g. IContesse et al. 
20041 : Iwilhelm et allbood) . instead of a decrease. Moreover, IPoUa & SolomonI (120081) reported detection 
of Alfven waves from the study of the variation of line widths in coronal hole off-limb region by analyzing 
SUMER/SOHO data, but they claim no evidence for damping of the Alfven waves because they found 
that the effect of stray light explained the observed decrease with height in the width of several spectral 
lines, starting about 0.1 -0.2 R© above the limb. This result was also confirmed from a later study by 
Dolla & SolomonI J2009h . 



The launch of th e Hinode mission opened new possibilities, thanks to the EIS spectrometer. Recently, 
Banerjee et al.l (120091) have identified the signature of Alfven waves in EUV spectra observed with Hin- 
ode/EIS spectrometer in a polar coronal hole through the increase of the nonthermal line-of-sight velocities 
up to - 150" above the limb. These waves appears to propagate undamped up to ~ 0.15 - 0.20 Rq above 
the limb, while it is at the moment unclear whether energy depo s ition i n the corona occurs above this al- 
titude. Similar work has been performed recently by iGupta et al.l (|2010l) . using spectroscopic observations 
of SOHO/SUMER and Hinode/EIS, where they conclude that the observed propagating disturbances are 
magneto-acoustic waves. Nevertheless, none of these works studied with EIS the behaviour of these waves 
at larger altitudes. 

Hence, in order to investigate the energy deposition by Alfven waves in polar regions above - 0.20 
Rq, between February 24-25, 2009 a special off-limb study of a polar coronal hole has been performed with 
the Hinode/EIS spectrometer (HOP 103, EIS study 340). As we describe here, these unique dataset allowed 
for the first time to study with EIS the line profile broadening up to - 0.40 Rq. The plan of the paper is 
as follows. In § 2, the observations acquired for this study and the data reduction techniques are outlined. 
In § 3, results of electron densities, non-thermal velocities and Alfven waves energy flux are presented. A 
discussion of the observational results and a comparison with theoretical studies are taken up in § 4, together 
with the final conclusions. 



2. Observations and data reduction 



2.1. EIS Study description 



The Hinode EUV Imaging Spectrometer (EIS; ICuUiane et all 120070 observes high-resolution spectra 
in two wavelength bands: 170-212 A (short wavelength band) and 246-292 A (long wavelength band). The 
instrument has 4 available entrance slits: two are available for imaging, with projected widths of 40" and 
266", and two are dedicated to spectroscopy with projected widths of 1" and 2", oriented along the North- 
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South direction. Pixel sizes on the detector are 1" (spatial dimension) and 0.0223 A (spectral d i mension). For 



single exposures the maximum field-of-view in Solar-Y is 512" (see iKorendyke et all 120061 : ICulhane et al 



20071. for a detailed description of the EIS instrument). 



The aim of the Hinode/EIS study # 340 was to measure the non-thermal velocities of various coronal 
ions above a polar coronal hole up to - 0.4 above the limb. Hence, between February 24-25, 2009, the 
512" long EIS slit was centered above the South polar coronal hole, 202" above the limb, covering a range 
of altitudes from - 0.05 R© below the limb up to - 0.48 R© above the limb (Fig. [Jl. In order to maximize 
the signal-to-noise ratio and perform reliable gaussian line fittings up to ~ 0.4 R©, the instrument acquired 
- 21.6 hours of continuous "sit & stare" observations with the 1" slit centered at X^O" and Y=-l 162" (155 
exposures, 500 s exposure time). These observations allowed us to derive the first EIS "deep field spectra" 
of EUV spectral lines in the off-limb polar coronal hole. Moreover, in order to have an idea of the structures 
(plumes or inter-plumes) present in the surrounding corona, two spatial rasters, each -1.19 hours long, have 
been acquired covering 80" in the X direction with the 2" slit (40 exposures, 100 s exposure time - hereafter 
"context study"). The first "context study" was performed before the "sit & stare" observations, the second 
one at the end of the "sit & stare" observations. The field of views (FOVs) covered with the spatial and 
temporal rasters are shown in Fig. [T] 

Due to Hinode telemetry limitations, for this study we were able to select only 10 spectral panels (32 
spectral bins per panel, equivalent to - 0.71A). These panels are centered over the following lines: O vi 
il84.00, Fe xu /1186.85, Ca xvu ^192.82, Ca xiv /1193.87, Fe xu il95.12, Fe xm ^202.04, He n .1256.32, 
Fe XIV /1264.78, Fe xiv A214.il, Fe xv /1284.16. Over the above spectral intervals, all the following Fe lines 
have been detected in this study with good statistics (i.e. with line peak intensity at least 2 times larger 
than the surrounding background): Fe x /1193.71, Fe xu /1186.88, Fe xu /1193.51, Fe xu A\95.\2 and Fe xm 
/1202.04. 



2.2. Data reduction and analysis 



Acquired data have been first calibrated with the standard routines provided by the SolarSoftware Hin- 
ode/EIS package, aimed also at removing instrumental e ffects such as the tilt of the EIS slits relative to the 
detector column's orientation, "hot" and "warm" pixels (|Youngll2009l. by interpolating missing data points; 
see), and the line centroid shifts related to the Hinode spacecraft Sun-syncronous orbit, occurring with a 
period of - 100 minutes. Second, by employing the period ic wavelength shift provided by the standard 
calibr ation routines (recently improved by the EIS Team; see lKamio et al.ll20ld : lKamio. Fredvik & Young 
201 ih . each line profile has been resampled bin per bin by fractional amounts of spectral pixels, preserving 
the total number of counts accumulated over the line profile (Fig.|2j left panels). Third, each identified spec- 
tral line has been fitted with a gaussian profile, in order to derive the evolution with time of the integrated 
emission at each altitude covered by the 5 12 slit spatial bins. Fourth, a further correction has been applied 
to the intensity versus time maps in order to remove a periodic pointing displacement by ± 6" observed in 
the Y direction, related to the spacecraft orbit and observed as an apparent periodic limb oscillation (Fig. |2] 
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right panels). 

Resulting intensities do not show clear transit of plumes/interplume structures dragged by solar rotation 
in and out of the spectrometer FOV (Fig. |3] right panel). This is in agreement with the global diffuse 
appearance of the off-limb EUV corona during our observations, as shown by the context EIS observations 
acquired before and after the "sit & stare" observations (Fig.|3l left panels), and also by FIT images (Fig.[l]l. 
No clear plume or inter-plume structures are visible at that time, and the line intensity evolution simply 
shows a very slow decrease during our observations. We were then unable to distinguish between these 
different structures. In order to have a better statistics, in this analysis we have summed over time aU data 
acquired during the "sit & stare" study. Hence, physical parameters we will derive in the next sections have 
to be considered as representative of an average plume/interplume plasma. Moreover, at higher altitudes, 
summing over - 21.6 hours of observations was not sufficient to have a good statistic on the line profiles. 
Hence, as also done by lBanerjee et al. ( 2009 ). profiles at the southern off-limb edge of the slit were averaged 
over different spatial intervals, summing over a larger number of pixels at higher altitudes. An example of 
line profiles obtained at 3 different altitudes with different spatial averages (considered in the analysis) is 
shown in Fig.|4l it is clear that the spectral lines are detected with a good statistics even at 0.4 R© above the 
limb. 



2.3. Stray light estimate 



The line profiles obtained in the Fe xii A195 spectral window (194.785-195.477 A) by averaging over 
the whole observation interval are shown in Fig. |5](top left panel) at 3 different altitudes. This plot shows 
that on-disk two more spectral lines are present, centered approximately at /i 194. 828 A and /1 195.41 6 A, as 
we derived with a multi-gaussian fitting. Both these lines disappear quickly moving off-limb, as shown by 
the average profiles at 50" and 150" off-limb in this Figure. The Fe xn /il95.119 is observed in our spectra 
at /1195.139, hence there is a shift by AA =0.020 A with respect to reference wavelenght. This means that 
the corrected centroids of the unidentified lines are /1194.808 and /1195.396. 

The identificati on of these two sp ectral lines is not straightforward: the latest version of the CHIANTI 
spectral code (v.7.0; IPere et all l2009h provides in this wavelenght interval only 4 observed spectral lines 
(Fe xu /1195. 119-^195.179, Ni xvi /1195.275 and Zn xx /1195.379) and 2 theoretical but unobserved hnes 
(Fe IX /1195.250 and Fe x /1195.316). Hence, the CHIANTI code does not prov ide any possible ide ntification 



Brown et al. 



(120081) identified 



for the /1194.808 and /il95.396 lines which we have observed. In the literature, 
with EIS in this range t he Fe xiv 1195.274 and Fe x /1195.395 lines (observed at the limb and 20" off-limb), 
while iLandi & YoungI (120091) identified a line observed at /1195.4 15 as the Fe yii /tl 95.391 line. Other 
possible lines close to the Fe xn /1195.12 have been reported by lYoung et al.l (l2009h (see their Fig. 6), 
which identi fied the Fe xiv /tl95.25, Fe x /1195.40 and Ni xv /tl95 .52; Ni xiv / 1195.2 7 was also reported by 
Young et alJ (|2009i ). A line observed at /1195.39 was identified by bel Zannal {200% as the Fe vm 1195.39 
line. Nevertheless, Fe xiv /1195.25, Ni xv 1195.52 and Fe vm 1195.39 were reported in Active Regions in 
these works. In particular, the Fe xiv and Ni xv lines have peak emission temperatures of ^ 2.0 x 10^ K and 
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2.5 X 10^ K, respectively, and no significant emission is expected from thiese lines in a coronal hole below 



0.4 Rq off-limb, \yhere SUMER and EIS data show electron temperatures tipically below 1.1 x 10^ K jLandi 



2008 



Hahn et al.ll201oi) . 



Hence, there are at least three possible candidates for the identification of the line we observed at 
/1195.396, i.e. Fe vm /1195.388, Fe x /1195.395 and Fe vii 195.391. The main difference between these lines 
is obviously the peak temperature of their emissivities: while the Fe x line emissivity peaks at T = 1.3 x 10^ 
K, hence at typical coronal temperatures, the Fe vu emissivity peaks at T = 4.0 X 10^ K (as provided by 
CHIANTI code), hence it is mainly a transition region line. This means that it is possible to distinguish 
between these two lines by looking at their intensity profiles with altitude. In fact, plots in the top-right 
panel of Figure |5]show that the intensities of Fe xii /1195.12, Fe xm /1202.04 and Fe x /1193.71 coronal lines 
have very similar decays off-limb, while the intensity of the unidentified line at /1195.396 decays much faster 
(in particular much faster than the Fe x /1193.71 line) and is mainly unobserved above - 120" off-limb (- 
0.13 Rq). This suggests us that the line we observe at /1195.396 is not a coronal line and could be identified 
as the Fe vu /1195.391 transition region line. Nevertheless, a contribution also from the Fe vm /1195.39 line 
cannot be totally excluded, in principle, because the contribution function of this line peaks at a temperature 
of - 6.3 X 10^ K and has a broader temperature distribution, making possible some contribution from the 
1.1 X 10^ K coronal hole plasma, even if (as far as we know), this line has been never reported off-limb in 
this kind of region. 

The identification of our unidentified line as Fe vn /1195.391 cannot be verified by comparing its inten- 
sity with those observed in other lines from the same multiplet, because no other Fe vn lines are present in 
our dataset. Given the observed Fe x /1193.71 intensity, we can at least estimate the expected contribution by 
the Fe x /1195.395 line to the observed unidentified line. As provided by the CHIANTI spectral code (v. 7.0), 
the Fe x /1195.39//1193.71 line ratio is almost constant in the off-limb coronal hole density range (10^ - 10^ 
cm"^) on the order of - 0.50. From the observed Fe x /1193.71 intensity we conclude that the expected Fe x 
/1195.39 intensity should be smaller than 12% and 43% of the total intensity of the unidentified line, 5" and 
390" above the limb, respectively. These values represent only upper limits, because have been computed 
before the correction for stray light, i.e. before removing the part of the radiation which is emitted on-disk 
and that, following optical paths not originally intended in the spectrometer, is finally detected off-limb. 
Hence we conclude that the Fe x /1195.39 coronal emission is expected to be much smaller than the intensity 
observed in the unidentified line. Unfortunately, we have no other Fe vm lines in our dataset, hence we were 
not able to perform the same estimate of the expected Fe vm /1195.39 contribution to the total intensity of 
our unidentified line. 

In summary, taking into account that, as mentioned above, a) the intensity of the unidentified line 
decays with altitude much faster than other coronal lines, b) that coronal emissions from other lines observed 
in this range (such as Fe xiv and Ni xv) are expected to be negligible in a coronal hole below 0.4 Rq off- 
limb, and c) that the Fe x /il95.39 coronal emission is expected to be much smaller than the intensity of the 
unidentified line, we conclude that the observed unidentified line mainly is not due to coronal emission, even 
if some contribution from Fe vm /1195.39 cannot be totally excluded. For these reasons, the line observed 
at /1195.396 has been identified here as the Fe vn /1195.391 transition region line. If this interpretation 
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is correct, its off-limb intensity should decay with altitude much faster than what we observed, becoming 
unobservable ~ 3 - 4" off-limb (corresponding to an altitude of ~ 2 - 3 • lO^km). This implies that the 
intensity we observe above this altitude can be ascribed to stray hght of the EIS instrument. 



As recently estimated by the EIS Team (|ugarte-Urra 2010 ). the stray light contribution for on-disk ob 



servations is on the order of 2% of the total intensity. I n order to es t imate the off-limb stray light contribution, 



we employed the same technique recently applied by iHahn et al.l (1201 ih : given the intensity lfevu{disk) of 



the Fe vii 195.39 line observed on-disk, we assumed that 2% of this intensity corresponds to the total Fe vu 
brightness scattered off in the bins looking in the off-limb corona. So, the stray light intensity is simply 
Ipevnistray) - 0.02 • Ifeviiidisk). The ratio between the Ifeviii^tray) and the Fe vu intensity observed 
off-limb at different altitudes h provides the stray light percentage curve htrayih), which is given by 



' stray 



(h) - 0.02 • lFeVIl{disk)lheVIl{h) (%). (1) 



Notice that the stray li ght computed with this technique may be slightly underestimated, because the 2% 



estimate made by Ugarte-Urral (120101) corresponds to th e "minimurn stray light level above the dark current" 



In any case, this estimate was also recently validated by lHahn et al.l(|2011h by checking the ratio between the 
Si X /1256. 3/261.0 lines observed immediately off-limb after removing from the He ii 1256.3 - Si vii /1256.3 
blend the 2% of the average He n disk intensity. The stray light curve resulting from the above equation is 
shown in the bottom right panel of Fig.|5](dotted line): the stray light contribution increases with increasing 
off-limb altitudes up to - 62 % at an altitude of 0.4 Rq off-limb. Even if, as explained above, the straylight 
values could be in principle underestimated, both because the 2% value is the minimum on-disk stray light 
level, and because a contribution from the Fe x /1195.39 and Fe vui /1195.39 coronal emissions cannot be 



totally excluded, this curve is in very good agreement with stray light values recently derived by iHahn et al 



(|201 ih at 0.10, 0.15 and 0.20 Rq with other spectral lines (triangles in the bottom right panel of Fig. |5]l. 
This suggests that the Fe x /il95.39 and Fe vm /1195.39 contributions to the Fe vu /1195.39 intrensity was 
negligible, at least up to 0.2 R©. 

The stray light curve Istray{h) as a function of altitude (Fig. [51 bottom right panel, dotted line) is well 
reproduced by the following function: 

htrayQi) = A arctau [B{h - C)] +mh + q (%) (2) 

with h measured in Rq from the limb (h > 0.03 Rq). By fitting the stray light curve with this function we 
obtained the following values for the above parameters: A = 0.214 (%), B ^ 19.5 (l/R©), C = 0.158 (R©), 
m = 0.153 (%/Ro), q = 0.269 (%). The resulting fitting curve (Fig. |5] bottom right panel, solid line) has 
been applied for the stray light correction of the intensity profiles of all the spectral lines. In particular, Fig. |5] 
(bottom left panel) shows the resulting off-limb stray light contribution (dashed line) to the observed Fe xu 
/1195. 12 intensity (solid line), together with the intensity profile corrected for stray light contribution (dotted 
line). Notice that the intensity profiles before the stray light subtraction (Fig. |5] top right panel) are quite 
flat above - 0.3 R©, suggesting possible stray light contamination, while after the stray light subtraction 
(Fig- El bottom left panel) the Fe xu /il95.12 intensity profile (dotted line) turns out to be slowly decaying 
even above that altitude. 
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3. Results 

The main target of this work is the estimate of the total energy flux F„ transported by Alfven waves 
off-limb in polar coronal hole up to - 0.4 Rq. To this end, we need both a determination of the electron 
density rig and the non-thermal velocity ^ at different altitudes. The determination of these three quantities 
(ng, ^, and F^) is decribed in the next Sections. 



3.1. Electron densities 



As mentioned above, over the available spectral intervals, we have detected and identified the Fe x 
.1193.71, Fe XII /1195.12, Fe xu /1193.51, Fe xii /il86.88, and Fe xm /1202.04 spectral lines. With these 
spectral lines, the better av ailable density dia gnostic technique is provided by the Fe xu /1195.12//1186.88 
density sensitive ratio (see I Young et al.ll2007i) . Electron densities obtained with this ratio are shown in 
Fig. [6] (left panel, open circles) up to an altitude of ^ 0. 12 Rp. These densitie s are in good agreement with 
those obtained from previous EIS off-limb observations (IBanerjee et al.ll2009l e.g.) for a polar coronal hole 
(Fig. [6l right panel, boxes and triangles). Nevertheless, above the altitude of - 0.12 R© the intensity of 
the Fe xu /II 86.88 starts to decay slower than the Fe xn /1195.12 intensity, resulting in a decrease of the 
/il95. 12//11 86.88 ratio, and finally an unrealistic increase of derived densities. The reason for this behaviour 
could reside in the blending of the Fe xu /1186.88 line (which is also a blend betwee n two Fe xu lines a t 
/il86.854 and .1186.887) with the nearby S xi /1186.84 hue, even if, as pointed out by lYoung et all J2OO9I) . 
the contribution of the S xi line to the Fe xu line should be negligible. 

In any case, for the above reason, we were unable to derive reliable electron density measurements 
above that altitude with the ratio /1195.12//1186.88. Unfortunately the other well detected spectral lines do 
not provide any other good density sensitive ratio. Hence, in this work we derived the electron densities 
above the altitude of 0.12 R© with another technique. For any spectral line, the CHIANTI spectral code 
allows the computation of the contribution function Gune(Te,ne,Ax) (erg cm^ s"^ sr~'), which depends on 
the unknown plasma electron temperature T^, density n^, and also on the unknown elemental abundance Ax- 
Given Gune, the observed intensity Inne of a spectral hue formed by only coUisional excitation can be written 
as 

hvUne C ^ rrr^ i \ . 2 j_ . -2 -1 



lline ~ 



J 

JLC 



Giine{Te,ne,Ax)nldz (erg cm ^s 'sr ') 



(3) 



47r J LOS 

where the integration is performed along the hne of sight (LOS). For spectral hues considered in this work 
the contribution from radia tive excitation has been neglected, as usually done in the so-called "coronal 
model approximation" (e.g. 



Mason et al. 



19971). The observed line intensity variations with altitude in the 



solar corona are due mainly to variations of n^ and T^, while the elem ental abundance Ax is usually assumed 
to be constant, in first approximation (see e.g. iDel Zanna et al. I I2OO2I) . Hence, given a pair of possible values 
for the average density and temperature along the LOS (< >, < Tg >), the expected line intensity 
can be estimated, by assuming a thickness along the LOS of the emitting plasma. By comparing the line 
intensity observed at a fixed altitude with all the expected values lune computed for a range of possible 
densities and temperatures, a curve on the (< rig >, < Tg >) plane of allowed values can be derived. Hence, 
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from a set of intensities for different spectral lines observed at the same altitude, it is possible to perform 
a minimum chi-square analysis and the best (< Ug >, < >) values have been derived by minimizing the 
difference between the observed and computed intensities for all the spectral lines (see e.g. Fig.|7]). 

In this work we have applied the above technique. In particular, for the computation we employed the 
intensities observed at different altitudes of the following iron spectral lines: Fe xu /1195.12, Fe xii il93.51, 
Fe xm /1202.04 and Fe xu /1 186.88. These lines have been selected as those observed with the better signal-to- 
noise ratio up to an altitude of - 0.4 R©. From the comparison in Fig.|6](left panel) it turns out that densities 
derived with our fitting technique are in quite good agreement with those derived from the ratio technique 
between - 0.03 and - 0.12 R©, while at lower altitude are probably underestimated. This difference close 
to the limb is likely due to our assumption that the component due to the resonant scattering is ne gligible, 
while it has been shown that this hypothesis is not always correct (see ^ Schrijver & McMuUenlboO O). Since 
the radiative and coUisional components are roughtly proportional to Ug and n^, respectively, neglecting the 
radiative component may lead to an underestimate of the density. On the contrary, densities derived w ith the 
Fe XII line ratio are not affected by the background radiation field (see discussion by lYoung et al ]l2009h . 



even 



if densities derived with thi s technique can be not accur ate close to the limb if the hypothesis of optically 



thin plasma falls down (see ISchrijver & McMullenll2000l) . Also, Fig.[6|( ^right panel) shows th at the densities 



values derived in this work are in good agreement with those derived by Banerje e et al.l(l2009h with EIS data 



(triangles and b oxes for plumes and inter-plumes, respectively), but slightly larger than densitie s derived by 



Wilhelml (120060 (upper and lower diamonds for plumes and inter-plumes, respectively) and by iDoyle et al. 



(Il999h (dash-dotted hne) with SUMER data. 



Uncertainties shown in Fig. [6] have been computed from the chi-square curves (Fig. |7J, starting from 
statistical uncertainties on the selected spectral line intensities. The latters are very small, because data have 
been averaged over - 21.61 hours of observation (i.e. smaller than 2% at each altitude). A power law fitting 
curve to the densities we derived above - 0.02 R© is given by 

n^h) = 5.644 • lO'^h'^^^-^ + 2.285 • lO^/z"^" '^^ + 5.656 • lO^/j"^'^" (cm"^) (4) 

with h (Rq) heliocentric distance measured from the Sun center. Because of the shape of the chi-square 
curves (see Fig.|7]), uncertainties in the resulting temperatures are larger than those on densities; hence we 
were not able to provide a reliable temperature profile at different altitudes. Resulting temperatures are on 
average T = (1.3 ± 0.3) • lO^K. This temperature is in agreement with the fact that for the above analysis 
we employed only line emissions from Fe^^"*" and Fe^^^ ions, whose maximum formation temperatures are 
around T 1.3 • lO^K. 

In order to verify the correctness of our results, the above densities and temperatures have been em- 
ployed to compute again the expected line intensities: a comparison between the observed (diamond sym- 
bols) and computed (solid l ines) intensities is shown in Fig. [8l In the above computation we assumed the 
Amaud & Roth enflug* (^1986) ionization equilibrium. Nevertheless, systematic uncertainties in the com- 



putation of expected intensities are mainly due to uncertainties in the atomic parameters and ionization 
equilibria assumed by the C HIANTI routine to compute the contribution functions, as also pointed out by 



Del Zanna & MasonI (|2005l ') for Fe xii lines. In particular, we estimate that different ionization equilibria 
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correspond to a systematic uncertainty by 24% and by 40% for the Fe xu and Fe xui lines, respectively. Cor- 
responding upper and lower limits for the computed intensities are shown as dotted lines in Fig.[8l this figure 
shows a general good agreement between the observed and computed line intensities at different altitudes, 
except for the Fe xii /^186.88 above - 0.12 R©, as already pointed out. 



3.2. Non-thermal velocities 

The non-thermal velocities of emitting ions have been determined from the observed line profiles as 
follows. The observed Une profile 1 le width, cr, for any coronal spectral Une can be written as 

where r,o„ is the kinetic temperature of the emitting ion with mass Mio„, ^ is the non-thermal broadening due 
to unresolved plasma velocities occurring along the line-of-sight and o"/ is the instrumental line broadening. 
This formula is valid under the hypothesis that ^ velocities have a Gaussian distribution. Given the observed 
line profile cr, the measured o"/ and by assuming a value for Tjnn, this relationship can be employed to 



estimate ^. As recently pointed out by the EIS Team ( Young 2011a)), the EIS instrumental width cr i is not 



constan t with different F-pixels along the slit. In this work we assumed o"/ values provided by lYoung 



(|2011al) . which are on the order of crj = 2.52 pixels (i.e. o"/ = 0.0563 A) and cr/ = 2.97 pixels (i.e. 
o"/ = 0.0663 A) respectively at the bottom and the top o f the EIS sUt window for our observations (Fig. [T). 
These values are slighly larger than the one derived by Brown et al. I JIoOB) (CT/ - 2.42 pixels - 0.054 A) 



for the short wavelength band by assuming ion kinetic temperatures r,o„ = 10^ K and non-thermal mass 
motions ^ = 30 km s~^ In any case notice that, given the uncertainties in the unknown quantities r,o„ and ^ 
(see later) and the relatively small variation in crj values along the slit given above, even in the case where 
CTj is over- or under-estimated, this will be mostly a systematic uncertainty, i.e. not affecting the relative 
variations of the derived ^ values at different altitudes. 

In this work we empl oyed the Fe xu /1195 .12 profiles, which give the better signal to noise ratio at all 



altitudes. As discussed by lYoung et all (120091) . the Fe xu A\95 line is a blend between two Fe xu lines at 
/1195.119 and /1195.179, resulting in general in asymmetric Fe xu A195 profiles when this line is observed 
on-disk. The line ratio Fe xu /1195.119//1195.179 depends on the electron densities, as shown in Fig. |9] 
(left panel). In particular, for coronal hole density values we derived off-limb above a polar coronal hole 
{ue ~ 10^ - 10^ cm"^; see Fig.|6]l, the ^195.179 Une is expected to be ~ 200 - 1000 times fainter than the 
/1195.119 line. Hence, the Fe xu A195 line is basically unaffected by the blend and the line profile can be 
safely fitted with a standard single gaussian profile off-limb above a polar coronal hole. 

Nevertheless, the same is not true for the average Fe xu /1 195 line emitted on-disk, which is clearly 
asymmetric (Fig. |9] right panel). In particular, for electron density values measured on-disk at the base of 



the corona on the order of 10*^-10" cm ^ (e.g. lYoung et al.ll2009l) the /1195.179 line is expected to be only 



~ 5 - 40 times fainter than the A195. 119 line. By fitting the average on-disk profile (Fig. |9j right panel, solid 
line) with two gaussians with the same width (i.e. the same Fe'^^ ion kinetic temperature) and separated by 
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0.06 A, it turns out that a /1195.119//1195.179 line intensity ratio by 5.91 (corresponding to a density on the 
order of - 4 x 10 " cm"^) well reproduces the asymmetry of the Fe xu A195 blend profile observed on- dislfl 

After a standard single gaussian fitting of the coronal line profiles, the cr values have been corrected 
for the instrumental profile broadening, as explained above. Corrected cr values provide an estimate for 
the broadening due at the same time to ion kinetic temperatures and plasma motions. In the literature, at 
least 2 different possible assumptions exist for the unknown quantity 7,0,,: the hypothesis of thermodynamic 
equilibrium {Tjo„ = Tg) and the hypothesis of ionization equilibrium (r,o„ = T,„ax)- Given the uncertainties 
mentioned above in the estimate of Tg, in this wo rk we consider the ion ization equilibrium hypothesis. By 



assuming (as also done by previous authors, i.e. iMazzotta et al.lll998l) that the unknown ion temperature 



equals the peak in the Fe xii /1195.12 line emissivity as provided by the CHIANTI spectral code (T = 
lO^'^K), we estimate the broadening, ^, due only to unresolved mass motions. 

Nevertheless, another effect needs to be taken into account. As mentioned above, the Fe xn /tl95.119- 
/1195.179 blend results in a broader Fe xn A195 observed profile emitted from the disk. This implies that, 
for this line, opposite to what typically occurs for other spectral lines, the stray light contribution results in a 
broadening of the coronal line profiles, and not in a narrowing. In order to correct the Fe xu /1195. 12 coronal 
profiles for the stray light broadening, we subtracted to the observed lines a broader stray light profile with 
a /1195.119//1195.179 line intensity ratio by 5.91 and relative total intensity given by the stray light curve 
derived in the previous Section. In order to quantify the effect of the stray light profile subtraction, line 
widths have been also derived with a single gaussian fitting before this subtraction. 

Resulting Fe xu line profile Full Widths at Half Maximum (FWHMs) and non thermal velocities ^ (km 
s~*) are shown in Fig.[TO](left and right panel, respectively), before and after the stray light subtraction (solid 
and dotted lines, respectively). Plotted FWHMs have been derived before the instrumental line broadening 
subtraction, while ^ velocities have been computed by removing both the instrumental and thermal broad- 
enings (formula 5). This figure shows that, within the error bars provided by the fitting procedure, a general 
trend is observed: after the stray light correction, the FWHM progressively increases up to an altitude of 
- 0.14 Rq above the limb, reaching a peak value of - 0.074 A, then decreases at higher altitudes down 
to - 0.061 A. Correspondingly, the estimated non-thermal velocity ^ increases up to - 40 km s~' at 0.14 
Rq above the limb, then decreases above that altitude down to - 10 km s~^ at 0.4 R© above the limb. The 
FWHM and ^ profiles shown in Fig. [10] are quite well reproduced by the following second order polynomial 
fitting curves: 

FWHM{h) = 7.075 • 10"^ 4.424 • 10"2/j- 1.720- 10"^/i^ (A) 

^(h) = 3.377 • 10-^ 8.638 • 10 /i - 3.565 • lO^/j^ (km/s) (6) 

(V) 

where h (Rq) is the altitude above the limb; the above fitting curves are shown by dashed lines in the left 
and right panels of Fig. [TOl 



also lYounj j2011bl) for a discussion of gaussian fitting and line blend for the EIS instrument. 
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3.3. Alfven waves energy flux 



In what follows we discuss the possible consequences of the working hypothesis that the unresolved 
plasma motions ^ we observed are due to the propagation of Alfven waves along open field lines above a 
polar coronal hole. Given the electron density and non-thermal velocity profiles derived above, it is possible 
to estimate the energy flux transfered into the corona by Alfven waves. In particular, the Alfven wave energy 
flux density f„ (e rg cm"^ s~^) for vertical propagation in a stratified, plane parallel atmosphere is given by 
(e.g. 



Moran 



200lh : 




(8) 



where v^-BI ^jAnp is the Alfven wave velocity, B is the field strength, and p is the mass density. This sim- 
plified relationship is correct only in the so called "Wentzel-Kramers-Brillouin (WKB) approximation", i.e. 
when the wavelengths are small compared to the background scale lengths. This approximation is valid on 
open fieldlines above ~ 0.03 - 0.04 R© off -limb, while non-WKB effects are important for wa ve propagation 
and reflection below this altitude (see e.g. lHollweglll981uCranmer & van Balle gooijenll2005l) . hence we can 
safely apply the above relationship for our data analysis. If the magnetic field flux 5 • A is conserved along 
opening flux tubes with cross sectional area A carrying the Alfven waves, and if Alfven waves propagate 
undamped, then the total Alfven wave energy flux F„ = f„A (erg s~^) should be a conserved quantity as the 
waves propagate outward. 

Values of f^, shown in Fig. [TT] have been computed at differen t altitudes with the above formula by 
assuming a magnetic field B = 8 G (as done by lBanerjee et al.ll2009l) . Results show a continuous decrease 
from f„ - 1.2 • 10^ erg cm"^ s"^ at 0.03 Rq down to f„ - 8.5 • 10^ erg cm"^ s"' at 0.4 Rq. In this altitude 
range a linear fitting to the /„■ decay provides an Alfven waves energy decay rate Af^/Ah - -4.5 • 10~^ erg 
cm"^ s~^ Below 0.03 Rq, basically because of the much stronger density gradient, decays much faster 
with a rate Af^jMi - -1.07 • 10~^ erg cm~^ s~^ Implications from these results are discussed in the next 
section. 



4. Discussion and Conclusions 

In this work we have analyzed special long duration sit & stare observations acquired by the Hinode/EIS 
spectrometer above the South polar coronal hole. The aim of this work is to study the EUV coronal line 
profile broadening observed off-limb and to provide a reliable measurement up to ~ 0.4 R© above the limb. 
In order to improve the signal-to-noise ratio, spectra acquired at different times have been averaged over 
the whole observation interval (~ 21.61 hours), and a spatial average has been also performed for spectra 
acquired at higher altitudes. Results show that the Fe xu A\95 line profile (the strongest we detected) 
progressively broadens going from the limb up to ~ 0.14 R© above it, reaching a peak FWHM value o f 



- 0.074A. This result is in agreement with what was previously obtained with EIS by lBanerjee et al.l(|2009h . 
as shown in Fig. [T2](diamonds). 

Thanks to the FOV of the EIS ofF-limb study we performed (Fig. [T), we were able also to study the 
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FWHM versus altitude profile up to further distances than those ever explored by previous authors with 
the same instrument. Our results show that above 0.14 R© the FWHM progressively decreases down to 
- 0.061 A at 0.4 R© above the limb. As we discussed in Section 2.3, the o bserved line width decrease a bove 
~ 0.14 Rq is not related to the EIS stray light contribution (as reported by lDolla & Solomonll2008L 12009. . for 
SUMER data), which has been estimated and removed from the observed Fe xii A195 line profiles. Hence, as 
a first result, we exclude that the line width decrease observed above ~ 0.14 Rq with EIS is an instrumental 
effect, and it needs to be physically explained. The main problem is that, as it is well known, the observed 
FWHM altitude profile is related to the altitude variation of two physical quantities (see formula 5), the ion 
kinetic temperatures r,o„ and non-thermal velocities ^ (after correction for the instrumental line broadening 
CT/ along the EIS slit), that can not be measured independently at the same time, unless some assumptions 
are made. Hence, derived non-thermal velocities are necessarily model-dependent. 

In the literature, many different assumptions have been made, hence providing in turn many different 
results. For instance: Baneriee et al. ( IQQsl) assu med a constant r,^,, = 1 • lO^K with altitude, in order to 



Tuetal 



(|1998r) assumed a constant ^ with altitude, different for different 

(|2003|) concluded that coronal 



Moran 



derive the ^ profile for Si ions; 

ions, in order to derive r,o„ as a function of the charge-to-mass ratios; 
ions have not a common velocity, hence non-thermal broadenings are not related to Alfven waves, but they 
assumed a common r,o„. In general, if a common temperature is assumed for different ions, there is no 
common non-thermal speed, and vice-versa. In this work we were able to study the line width variation with 
altitude up to ~ 0.4 R© only for the Fe xu A195 line, while other spectral line profiles have a much lower 
signal to noise ratio. Non-thermal velocity profiles have been derived by assuming a constant r,o„ = 10^-^ 
at different altitudes, i.e. equal to the Fe xii /1195 line maximum emissivity, as provided by the CHIANTI 
spectral code. After a correction for the instrumental line broadening, this temperature provides us a pro- 
gressive non-thermal velocity ^ increase with altitude above the limb up to - 40 km s~^ at 0.14 R©, then a 
decrease above that altitude down to - 10 km s~^ at 0.4 R©. 

In this work we also derived a coronal hole electron density profile for the off-limb corona up to 
0.4 Rq (eq. 4) with two methods: the line ratio technique and a minimum chi-square method aimed at 
reproducing the observed line intensity decays with altitude. The density and non-thermal velocity profiles 
derived led us to conclude that the Alfven wave energy flux density f„ continuously decreases with altitude 
above the limb from /„ ^ 1.2 • 10^ erg cm~^ s~^ at 0.03 R© down to /„■ - 8.5 • 10^ erg cm~^ s^^ at 
0.4 Rq, with an average decay rate by -4.5 • 10"^ erg cm"^ s"^. Hence, the Alfven waves energy flux 
density is not conserved (Fig. fTTl) . not even in the lower coronal hole region where the line profile FWHM 
increases with altitude (Fig.fTOl). The disagreement with the behaviour expected for undamped Alfven waves 
becomes more evident at larger altitudes. In fact, in the hypothesis of undamped Alfven waves propagation, 
densit ie s and non-ther mal velocities should be related by ^ oc p"'/^, as already pointed out by iHoUwe^ 
(119901) ; iMoranI (1200 lb . Hence, given a density profile p{h), it is possible to compute the corresponding 
non-thermal velocity profile ^(h) that should be observed for undamped w aves. Three examp le ^(h) curves 
are sh own in Fig. [T2lbv as suming coronal hole density profiles provided by Doyle et al.l (|l999h (dash-dotted 
line), Guhathakurta et al. ([l999) (dotted line) and the density profile derived in this work (dashed line). 
A comparison between these curves and non-thermal velocities for Fe^^"*" (diamonds), Fe^^"*" (triangles) 
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and Fe'^""" (circles) ions measured with dilferent techniques in a coronal hole show that, even taking into 
account the broad scattering of different measurements (hence the large uncertainties), non-thermal velocity 
measurements could be in agreement with the propagation of undamped Alfven waves up to no more than 
~ 0.2 Rq, while at higher altitudes the disagreement with this hypothesis is evident. 

The total Alf yen wave energy flux repo rted here is larger than what required to heat the quiet corona (3 • 



Withbroe & Noveslll977l). Ney ertheless. a physical explanation for this wave dissipation 



10^ erg cm"^ s~^; 

is needed. For instance, IZaqarashvili et al.l (l2006l) demonstrated that Alfven and sound waves are decoupled 
for values of the plasma yS parameter < 1 , but become coupled when /? = 1 , allowing energy transfer from 
Alfven to sound waves. If this energy is transferred to sound waves above the chromopsphere, then these 
waves will quickly dissipate their energy in the low density corona by inducing shock w aves. This c ould 
be a possible explanation for the observed Alfven wave damping: for instance following iGaryl (|200l b the 
plasma p above an active region approaches unity around 0.2 Rq, while outsid e active regions can approach 
unity even at lower heights. Alternative explanations have also been proposed. iDwivedi & Srivasta va (2006) 
showed that, by taking into account the effect of viscosity and magnetic diffusivity in the radial propagation 
of Alfven waves, an Alfven wave energy flux density decrease is expected in polar coronal holes above 0.21 
R©; this decrease is more pronounc ed for low frequency (10^ H z) compared to high frequency (10^ Hz) 
waves. This effect was confirmed by ISrivastava & Dwivedil (l2007h also for the rnore g eneral case of oblique 
propagation. We notice here that the decay obtained by iDwivedi & Srivastaval (12006) for the energy flux 
density of Alfven waves is in qualitative agreement with our results (Fi g. [TTIl in the case of dissipat ion due 
to magnetic diffusivity alone for low frequency (< 10^ Hz) waves (see IDwivedi & Srivastaval l2006l Fig. 2, 
middle panel). 

More recent observational results demonstrate the existence of ubiquitous low frequency (10~^ -2- 10"^ 
Hz) Alfven waves propagating th roughout the quiescen t atmosphere with sufficient energy to accelerate the 
solar wind and heat the corona dMcIntosh et al.ll201lh and suggest that low-freque ncy Alfvenic motions 
may suffer from significant damping as they propagate along magnetic structures dTomczyk & Mcintosh 
20091) . Nevertheless, other theoretical wor ks suggest that Alfven waves should p ropagate undamped at these 
low altitudes above the limb. For instance, ICranmer & van BallegooijenI (|2005l) developed a comprehensive 
model for the propagation of Alfven waves in a coronal hole from the photosphere to the interplanetary 
medium, and demonstrate that the frequency-integrated velocity amplitude should increase monotonically 
with height from the chromosphere up to 10 Rq. In the present work we demonstrate that the FWHM decay 
of EUV lines (reported also by previous authors) cannot be simply ascribed to instrumental effects (i.e. stray 
light). In principle, alternative physical explanations for this observational effect, other than the damping 
of Alfven waves, could be proposed. But a challenge remains to understand how to explain on the one 
hand the non-thermal velocity decay of heavy ions reported here in the EUV low coronal emission {h < 0.4 
Rq), and on the other hand the very large non- thermal velocities of heavy i ons discovered in the UV spectra 
emitted by the intermediate and extended corona {h > 0.5 R©; see review by Kohl et aPboOfi , and references 
therein). 
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Fig. 1. — Left: the EUV corona as seen by the SOHO/EIT telescope on February 24, 2009 with the Fe xii 
/1 195 filter. The white box shows the Hinode/EIS spatial raster's FOV, while the yellow vertical line shows 
the slit position during the "sit & stare" study. Right: a zoom on the region surrounded in the left panel by a 
yellow box. The EIS slit, extended off-limb out to ~ 0.48 R©, is shown as yellow line and the spatial raster's 
FOV is illustrated by the white box. 
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Fig. 2. — Left: Fe xii /1195.12 line profiles (x-axis) averaged along the spatial direction (512") as a function 
of time (j-axis) before (left) and after (right) the spectral resampling performed in order to remove the 
wavelength shift related to the Hinode spacecraft orbit. Right: a zoom over the Fe xu A195. 12 line intensities 
observed around the limb (j-axis) as a function of time (;c-axis) before (left) and after (right) the spatial 
resampling performed in order to remove the pointing shift related to the Hinode spacecraft orbit. 
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Fig. 3. — The Fe xii /il95.12 line intensity as observed during the context and the sit & stare observations. 
Distances are measured from the solar limb. Left panels (a-b): the line intensity as observed with two spatial 
rasters each ~ 1 hour long, covering 80" in the X direction with the 2" slit, acquired before (a) and after (b) 
the sit & stare observation. Right panel (c): the line intensity evolution during the ~ 21.6 hours of "sit & 
stare" observations with the 1" slit. 
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Fig. 4. — Example of line profiles for the Fe xii /1195.12 (top) and Fe xm /1202.04 (bottom) spectral lines 
at 3 difi'erent altitudes above the limb, as indicated at the top of each plot. Line profiles shown here were 
obtained by averaging line profiles over the whole sit & stare observation interval and by further averaging 
over 125 (right column) and 25 (middle column) arcsecs along the slit. 
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Fig. 5. — Top left: the spectral profiles observed in the Fe xii /1195.12 panel 50" on-disk (solid line), 50" 
off-limb (dashed line) and 150" off-limb (dotted line) from the solar limb. The solid line profile shows 
the Fe xu /il95.12 spectral line and 2 other lines, centered approximately at /il95.41 and /il94.83 (see text 
for discussion on the identification of these lines). Top right: normalized intensity profiles with altitude 
(logarithmic scale) for the Fe xu /1195.12 (solid), Fe xm /1202.04 (dotted), Fe x /1193.71 (dash-dots) lines 
and for the unidentified line at /il95.41 (dashed). Bottom left: intensity profile with altitude (logarithmic 
scale) for the observed Fe xii /1 195. 12 line (solid), the fraction of Fe xu emission from the disk (dashed), 
and the Fe xn intensity corrected for stray light (dotted). Bottom right: the EIS stray light fraction estimated 
with the int ensity of the /t 195 .41 line, identified as Fe vu /1195.415 (dotted line), compared with the values 
obtained by lHahn et al.l (1201 ll) (triangles). Solid line shows the fitting curve (eq. 2) of the derived values. 
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Fig. 6. — Left: comparison between electron densities derived from the emission measure fitting procedure 
(filled circles) and those derived from the Fe xu /1195.12//1186.88 line ratio (open circles). Dotted line 
shows the power law fitting curve (eq. 4). Right: electron densities derived with the emission measure 
fitting procedure ( filled circles) and the power l aw fitting curve (dotted line). Density profiles derived by 
Doyle et all (11999 ) (dash-dots), WiUiehnl (120060 (upper and lower diamonds for plumes and inter-plumes, 



respectively) and lfianerjee et all (|2009h (triangles and boxes for plumes and inter-plumes, respectively) are 
also shown for comparison. 
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Fig. 7. — Example of the chi-square distribution over the electron density vs. electron temperature plane for 
the line intensities observed at an altitude of 0.1 Rq above the limb; solution point at that altitude is shown 
by the white cross. 




Fig. 8. — Comparison between the observed line intensities (diamonds) and those computed (solid line) with 
densities and temperatures derived from the emission measure fitting procedure; upper and lower limits are 
shown as dotted Unes (see text). 
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Fig. 9. — Left: ratio between the emission measures (y-axis) of the Fe xii /1195.119 and /1195.179 lines 
provided by the CHIANTI spectral code (v.7.0) as a function of densities (x-axis). Right: Fe xii /{195.119- 
/1195.179 blend profile observed on-disk (solid), the corresponding two gaussians fitting curve (dashed), and 
single gaussians reproducing each spectral line (dotted). 
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Fig. 10. — Left: altitude profile of Fe xii /1195.12 FWHMs without correction for instrumental broadening 
as derived with a single gaussian fitting before (solid line) and after (dotted line) subtraction of the Fe xii 
stray light profile. Right: corresponding non-thermal velocities ^ corrected for instrumental broadening 
computed before (solid line) and after (dotted line) subtraction of the Fe xu stray light profile. Dashed lines 
show parabolic fittings to the dotted curves (see text). 
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Fig. 11. — The Alfven wave energy flux /„, as a function of altitude, computed by assuming magnetic flux 
conservation with non-thermal velocities ^ derived before (solid line) and after (dotted line) subtraction of 
the stray light profile. This plot shows that the Alfven wave energy flux progressively decays with altitude, 
with an average decay rate by -4.5 • 10"^ erg cm"-' s"' above 0.03 Rq. 
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Fig. 12. — Comparison between non-themial velocities ^ derived in this study (solid line) and those derived 
by several authors with different instruments above a polar coronal hole for the following Fe ions (in the 
on-line color version same colors corresponds to same ions): Fe'^"*" deriv ed with Hino de/EIS (diamonds; 



Banerjeeet all 120091) . Fe''+ derived with SOHO/SUMER (square boxes; 



Moran 



2003h. F e^^+ (triangles) 



and Fe'^"*" (circles) derived with spectroscopic total eclipse observations dSingh et al.ll201lh . These values 
are compared with tho se expected for prop agation of undamped Alf ven waves by assuming density profiles 
by bovle et all Jl999l) (dash-dotted line), louhathakurta et all (fl999i) (dotted line) and the density profile 
derived in this work (dashed line). 



